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We report measurements of absolute transition frequencies and hyperfine coupling constants for 
the 8S1/2 , 9Si/2 , 7D-J/2 , and 7D5/2 states in ^"^"^Cs vapor. The stepwise excitation through either the 
6Pi/2or 6P3/2 intermediate state is performed directly with broadband laser light from a stabilized 
femtosecond laser optical-frequency comb. The laser beam is split, counter-propagated and focused 
into a room-temperature Cs vapor cell. The repetition rate of the frequency comb is scanned and 
we detect the fluorescence on the 7Pi/2,3/2 — >■ 6S1/2 branches of the decay of the excited states. The 
excitations to the different states are isolated by the introduction of narrow-bandwidth interference 
filters in the laser beam paths. Using a nonlinear least-squares method we find measurements 
of transition frequencies and hyperfine coupling constants that are in agreement with other recent 
measurements for the 8S state and provide improvement by two orders of magnitude over previously 
published results for the 9S and 7D states. 

PACS numbers: 42.62.Fi, 42. 62. Eh, 32.30.-r, 32.10.Fn 
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I. INTRODUCTION 



The basic ideas of using repetitive pulses and mode- 
locked lasers for high-resolution precision spectroscopy 
were originally developed in the 1970's, with particular 
focus on Doppler-free two-photon transitions [iHaL- The 
advent of fully stabilized optical frequency combs 's'l in 
2000 renewed interest in many spectroscopic applications 
of mode- locked lasers (for a review see Ref . [9| ) . In many 
instances the use of direct frequency-comb spectroscopy 
(DECS) offers many advantages over spectroscopy done 
with cw lasers. In particular, the spectral versatility pro- 
vided by modern frequency combs allows the simultane- 
ous study of multiple atomic transitions with a single 
experimental setup. Additionally, DECS allows for a rel- 
atively simple determination of the absolute frequency 
by providing a direct link to the S.I. second. In this work 
we exploit this versatility to study multiple transitions in 
cesium (Cs). 

Precision spectroscopy of Cs plays an important role 
in atomic and fundamental physics, including measure- 
ments of atomic parity violation [13-113 1 searches for the 
permanent electric dipole moment of an electron [ij, [ij] , 
photon recoil measurements that help determine the fine 
structure constant [H, [l3| , measurements of the distri- 



bution of nuclear charge and magnetization [17| , and nu- 
merous ap plic ations in laser cooling, Bose-Einstein con- 
densation |l8|. atomic clocks [l^, and magnetometers 
[23. 

The present work focuses on measuring multiple tran- 
sition frequencies and the hyperfine structure of select 
low-lying states in cesium via two-photon excitation. 
The relatively simple experimental setup uses a room- 
temperature Cs vapor cell and a self-referenced optical 
frequency comb to excite Doppler-free two photon tran- 
sitions from the Cs ground state to the 8S1/2 , 9Si/2 , 
7D3/2 , and 7D5/2 states. Our results are consistent with 
previous measurements for the well-studied 8S1/2 state 
|2ll [2^ and provide significantly reduced uncertainties 
for the 9Si/2 , 7D3/2 , and 7D5/2 states 

This paper is structured as follows: Section llll contains 
a description of our experimental system and measure- 
ments. In Section [IIII we examine the features of the fiu- 
orescence spectra and discuss the theoretical model used 
to describe the results. Section IIVI describes the analy- 
sis of the two-photon transitions and the recipe used in 
extracting the transition energies. This section also in- 
cludes our results and the details of the systematic effects. 
We conclude our work in Section Ivl 



II. EXPERIMENTAL SETUP 
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The relevant Cs energy levels are shown in Eig. [T] 
The excitation pathways to the four states 8S1/2 , 9Si/2 , 
7D3/2 , or 7D5/2 proceed through a stepwise excitation 
through either the 6Pi/2 0r 6P3/2 states. We detect the 
two-photon excitation as blue fluorescence emission as 
the atoms decay from the 7Pi/2,3/2 states to the ground 
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FIG. 1: (Color online) Energy levels of "^Cs (/ = 7/2) that 
are relevant to this work. The excitation pathways are shown 
in solid red along with the transition wavelengths. The decay 
channels used in the detection are shown in blue dashed lines. 
The total angular momentum of the hyperfine components is 
listed next to each state. 



state. 

The experimental set-up is shown in Fig. [2] We ex- 
cite the cascaded two-photon transitions using a mode- 
locked Ti:Sapphire laser that generates a broad optical 
frequency comb. The portion of the output spectrum 
that was used extended from 650 nm to 1000 nm, cov- 
ering all of the excitation transitions of interest in Fig. 
[TJ The laser is described in detail in Ref. and here 
we only briefly mention the information relevant for the 
experiment. The Kerr-lens mode-locked Ti:Sapphire ring 
laser produces an output spectrum consisting of discrete 
equally spaced coherent optical modes. The frequency 
of each optical mode, Vn, can be related to two mi- 
crowave frequencies, the laser repetition rate, frep, and 
the carrier-envelope offset frequency, /o, by 



1^71 — n frep + fo 



(1) 



where n is an integer mode number 4 x 10^. The 
carrier-envelope-ofFset frequency is stabilized by use of an 
/-to-2/ interferometer, thus "self-referencing" all modes 
of the comb f?\. For the measurements described here 
fo was stabilized to 45 MHz. The repetition rate of the 
laser is « 1 GHz. The repetition rate was phase-locked 
to a signal generator that was referenced to a hydrogen 
maser. With the comb stabilized in this way, the frac- 
tional instability of every optical frequency of the comb 
is w 2 X 10~^'^T~2^ where t is the averaging period in 
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FIG. 2: Block diagram of the experimental set up. The com- 
ponents are: BSl, BS2 nonpolarizing beam splitters; Fl, F2, 
F3 interference filters; L lens; PD photodiode: PMT photo- 
multiplier tube; VC: Cs vapor cell. 



seconds. 

The laser beam was split by use of a 50-50 nonpolar- 
izing beam-splitting cube (BS 2). The two output laser 
beams were linearly polarized along the same axis. The 
beams were filtered by use of either narrow-band inter- 
ference filters or a combination of long-pass and short- 
pass filters to select specific transition pathways. The 
counter-propagating beams were then focused by lenses 
of 25 cm focal length to a spot size of ~ 100 /Ltm inside 
the vapor cell. Focusing the beams to a smaller waist size 
increases the two-photon excitation rate, but also leads 
to broadening of the resonances due to the finite transit- 
time of the atoms through the beams. The transit-time 
broadening for a waist of this size is about 2 MHz for 
room-temperature Cs atoms. The average optical inten- 
sities in the vapor cell were tens of W/cm^ spread over 
« 10 nm bandwidth, corresponding to « 4000 optical 
frequencies. 

An evacuated Pyrex Cs vapor cell was used in the ex- 
periment. The cell has a diameter of 2.5 cm and is a 
wedged shape with an optical path length ranging from 
0.5 cm - 1 cm. Windows were fused directly to the cell at 
an angle. The cell was carefully cleaned and filled on an 
oil-free vacuum system and was evacuated to a base pres- 
sure of ~ 10~^ Pa. A small amount of Cs was introduced 
into the cell prior to removing it from the vacuum system. 
The cell temperature was varied by resistive heater wires 
wound around the cell. We note that we observed shifts 
in the transition frequencies as large as 0.5 MHz in a Cs 
vapor cell that was not cleaned and vacuum processed as 
carefully as the cell used in these measurements. 

Most of the measurements in this experiment were 
taken with the Cs cell temperature at 297 K, with fluctu- 
ations less than 1 K. Data were taken at higher temper- 
atures to evaluate systematic shifts arising from pressure 
shifts as described in Section llVBI 

A set of six coils was used to zero out the earth's mag- 



netic field in the region of tlie vapor cell. The field a 
reduced to a level of < 1 /iT over the interaction reg 
inside the vapor cell. 

Two-photon resonance signals were monitored by m 
suring the 7Pj -> 6S1/2 (J = 1/2, 3/2) cascade fl 
rescence at 459 nm and 455 nm. The fluorescence a 
collected at a right angle to the two counter-propagat 
beams and was detected by a photomultiplier tube. T 
detector was shielded from background light and the sc 
tered light from the excitation laser beams by a 45 
nm filter with a 10 nm pass band (F3). The detec 
output was amplified and filtered (low-pass at 1 kHs 
dB/octave) before being sampled and recorded wit] 
computer-controlled data acquisition board (DAQ). 

The DAQ used in these measurements is similar 
that used in Refs. [H, [l^llll. The computer control 
the scan of the repetition rate by stepping the freque: 
of the signal generator to which the repetition rate i 
phase locked. Data scans were collected by changing 
frequency of the signal generator from 1000.771 725 M 
to 1000.774 725 MHz in 1 Hz steps. This 3 kHz rai 
in the repetition rate results in a 2-3 GHz variation 
the two-photon frequency, depending on the transiti 
At each 1 Hz step the fluorescence was digitized w 
multiple samples and averaged for 1 s. The statist: 
mean and the uncertainty in the mean were evaluatec 
each repetition rate frequency. Spectra were collected 
both increasing and decreasing frequency scans, result 
in 6002 data points. 



III. OBSERVED SPECTRA AND 
THEORETICAL CONSIDERATIONS 

Figure [3] shows the spectrum collected as the repet 
tion rate of the frequency comb was scanned withov 
using the filters marked as Fl and F2 in Fig. [2l Th 
spectrum consists of a multitude of narrow Doppler-fre 
peaks and a Doppler-broadened background. The narro' 
peaks are due to excitation to the 8S, 9S and 7D state 
by two counter-propagating photons, while the Dopple: 
broadened background is due to excitation from two cc 
propagating photons. The observed spectrum is quit 
complicated because of the numerous discrete frequencic 
present in the comb combined with the multiple atomi 
energy levels and transition pathways accessible. 

The spectrum is greatly simplified and the Dopple 
background is eliminated by introducing filters Fl an 
F2, as shown in Fig. [5] The presence of the filters r( 
stricts the transitions and enables the selection of ind: 
vidual excitation pathways. Figure |4] shows the spec 
trum observed when one of the counter-propagating laser 
beams is restricted to be between 845 nm to 855 nm and 
the other between 780nm to 800 nm. With this combi- 
nation, the first beam can excite the atoms on only the 
6S1/2 6P3/2 transition at 852 nm (the D2 transition) 
and the second beam can only excite the atoms on the 
6P3/2 ~^ 8S1/2 transition at 795 nm. This eliminates ex- 
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FIG. 4: The calculated spectrum (upper trace) and exper- 
imental signal (lower trace) of the 7P — > 6S fluorescence 
when the Cs atoms are excited to the 8S1/2 state through the 
6P3/2 state. The groups of peaks in the calculated spectrum 
are labeled by F-F" to denote the hyperfine ground state and 
hyperflne flnal state of the transition. The peaks in each group 
correspond to the different hyperfine components of the inter- 
mediate state. The peaks in the experimental spectrum that 
are labeled with letters were used in extracting the transition 
frequencies, as described in the text. 



citations to the 9Si/2 , 7D3/2 , and 7D5/2 states, as well as 
the excitation to the 8S1/2 state through the 6P1/2 state 
and removes the Doppler-broadened background. The re- 
sulting spectrum consists of groups of two or three peaks. 
The peaks labeled d, e, and / form one such spectral 
group and are shown in more detail in Fig. [51 Each spec- 
tral group is the result of excitation from a given ground- 
state hyperfine level to a given excited-state hyperfine 
level, and the individual peaks in the group correspond 
to specific transition pathways through every allowable 
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hyperfine level of the intermediate state. For example, 
the d, e, and / peaks in Fig. [5] all correspond to ex- 
citation from the the 6S1/2 {F = 4) ground state to the 
8S1/2 {F" = 4) excited state. The three peaks correspond 
to excitation through the 6P3/2 {F' — 3, 4, 5) intermedi- 
ate states. We see the groups recur as a different pair of 
optical modes is resonant with a given transition. This 
repeating pattern occurs at intervals of (5 /rep ~ 1-2 kHz, 
corresponding to a change in the two-photon optical fre- 
quency of w 1 GHz (i.e., the laser mode spacing). For a 
given group of three peaks there are three different veloc- 
ity classes that satisfy the resonance condition for each 
of the corresponding hyperfine components of the inter- 
mediate level. For another group of three peaks, the res- 
onant condition for the components of the intermediate 
level is satisfied by a different pair of comb modes, and 
by different velocity classes. Thus, the amplitudes of the 



peaks arising from the same ground, excited, and inter- 
mediate states are not the same in the different spectral 
groups. 

We can understand and calculate the spectra using 
the standard two-photon transition probabilities deriv- 
able from second-order time-dependent perturbation the- 
ory. We start by considering a transition of an atom 
in the ground state with total angular momentum F^ 
|6Si/2 F'), to an excited state with total angular momen- 
tum F'\ \n" L'jn F") via excitation by two cw lasers with 
angular frequencies wi and W2 , intensities Ii and I2 , and 
wave vectors ki and ^2- Here n" refers to the principal 
quantum number of the excited state. We further as- 
sume that (jji is close to the resonant frequency for the 
6S1/2 6Pj' transition, with J' = 1/2 or 3/2. The 
probability of the atom to make a transition is 26] 



P{6S,/2F, n"L'}„F") <x 



1 



2F + 1 



F:n" L'',F" 



OJl + kl ■ V 



^2 / 7„ 
LJ2 + k2 ■ v)\ + 



E 



E 

F'M' 



{n''L';,nF''M'^\e2-d\6Pj,F'M'p){'oP,j,F'M'p\ei ■ JjeSi/aFAfF) 



W6Si/2-F:ePj,F' - [l^l + fci ■ u 



76P J, 



(2) 



where v is the velocity of the atom, Mf, M'p, and Mp are 
the projections of the total angular momenta F , F' , and 
F" along the axis of quantization, d is the electric dipole 
operator, 61(2) is the unit vector along the direction of 
the polarization for the first (second) laser beam, "fnLj 
is the homogeneous line width of the state \nLj), and 
'^nLjF-.n' L'j,F' IS thc resonant angular frequency of the 
transition \n Lj F) — > \n' L'j, F'). For a given light polar- 



ization, specified by the vector component in the spheri- 
cal basis g, the matrix elements for thc different magnetic 
sublevels can be related to reduced matrix elements that 
do not depend on the magnetic sublevels or the total 
angular momentum (F, F', or F") by use of the Wigner- 
Eckart theorem and standard Clebsch-Gordan relations 
as \21\ 



{n'L'j,F'M'p\dq\nLjFMF) = (-1 



,F'-M' 



F' 1 F 

~M'p q Mp 



{n'L'j,F'\\d\\nLjF) 



= {-lf'-^''^+'+^+-J+^^{2F + 1){2F' + 1) 



F' 1 F 

-M'p q Mp 



J' F' I 
F J 1 



{n'Lj,\\d\\nLj). 



(3) 



Here I — ^ is the nuclear spin and the ( ) and 

I } terms are the usual 3- J and 6- J symbols. Thus, 
the transition probabilities for different ground-state an- 
gular momenta F and excited-state angular momenta F" 
can be related to a common factor by use of the appro- 
priate Clebsch-Gordan coefficients. If we assume the two 



laser beams are counter-propagating then examination of 
Eq. [5] reveals that for a given transition from the ground 
state with angular momentum F, through an interme- 
diate state with angular momentum f, to a final state 
with F" the transition probability will be a maximum 
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when 

Wl(l + ^) + ^2(1 - ^) ^ (^6Si/2F:7i" L'^„F" (4) 

and 

+ = W6Si/2-F:6Pj,-F'- (5) 

If we now consider the effect of having a number of dis- 
crete frequencies, coming from the different modes of the 
comb, we see that, as a result of the detuning from the 
intermediate state, only two near-resonant comb modes 
will contribute substantially to the transition probabil- 
ity. Because the spacing of the frequencies for the comb, 
frep ~ 1 GHz, is much larger than the homogeneous 
linewidths of the states, 5^ ~ 5 MHz, we are justified 
in considering only the two modes closest to the reso- 
nance conditions of Eqs. ^ and ([5]). Since these twc 
frequencies are coming from the comb, they are relatec 
to the offset frequency and repetition rate of the comf 
by two integers as 

ll>i =27r(/o + ni frep) (6; 

LO2 = 27r(/o + n2 frep)- {7] 

Substituting these relations into Eqs. and (O wf 
see that for a given pair of modes, ni and n2, and £ 
fixed value of the offset frequency there exists a unique 
repetition-rate frequency and velocity that will simul- 
taneously satisfy both resonance conditions. Therefore 
given an atomic sample with a broad enough velocity dis- 
tribution, there will exist a repetition-rate frequency thai 
will result in a stepwise resonant excitation of the two- 
photon transition for a specific velocity class. This leads 
to the narrow Doppler-free peaks seen in Fig. 21 Because 
the velocity class excited and the resonant repetition-rate 
frequency both depend on the energy of the intermediate 
state, the transition resonance for each different interme- 
diate state occurs at different repetition rates. Because a 
different velocity class contributes to the different peaks, 
the amplitudes of the peaks vary greatly, depending on 
the number of atoms that are present in the resonant 
velocity class. 



IV. ANALYSIS AND RESULTS 

A. Analysis 

We calculated the excitation probabilities by numer- 
ically integrating Eq. [5] over a Gaussian velocity distri- 
bution for optical modes near the resonant condition. A 
total of ~ 15 optical modes were used (10 for the first 
stage excitation and 5 for the second stage excitation). 
The calculated spectra agree quite well with the experi- 
mental spectra (Fig. 2]). Slight variations in the relative 
amplitudes of the peaks are attributed to slight deviation 
of the polarization of the excitation light from linear. 

The resonant frequencies for a transition from a ground 
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FIG. 5: The 7P — >■ 6S fluorescence when the Cs atoms are 
excited to the 8S1/2 state through the 6P3/2 state. The peaks 
correspond to the d, e, and / peaks in Fig. [l] This data set 
was taken at a temperature of 319 K and was used in evalu- 
ating systematic efi'ects related to temperature dependence. 



related to the center-of-gravity frequencies and the hy- 
perfine coupling constants through (see e.g. {2^ 1 



\k{F", J" f + Ik{F", J") - I{I + l)J"{J" + 1) 



2/(2/ - 1) J"(2J" - 1) 



(8) 



where k{F, J) = F[F + 1) - /(/ + 1) - J( J + 1). 

The repetition rates at which the calculated resonances 
occur depend on the center-of-gravity frequencies for the 
two stages of the excitation and the hyperfine A and B 
constants of the ground, intermediate, and final states. 



Spectra were calculated with different values of the two- 
photon center-of-gravity frequency and hyperfine cou- 
pling constants of the final excited state. The interme- 
diate state energy and coupling constants used in the 
calculated spectra were fixed to the values determined 
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from Gerginov et al. [l^ H^. The two-photon tran- 
sition center-of-gravity frequencies f6Si/2:n"L'^„ and the 
hyperfine coupHng constants ^„"l" and Bn"L" were 
extracted from the experimental data by least-squares 
comparison of the experimental resonant repetition-rate 
frequencies with the resonant repetition-rate frequencies 
determined from the calculated spectra. 

The resonant repetition rate for each of the peaks in 
both the experimental and calculated spectra were deter- 
mined from a least-squares fit. Fitting the spectra was 
necessary in order to account for the overlap of the dif- 
ferent peaks. The calculated spectra were assumed to 
have a Lorentzian line shape. However, the experimental 
spectra were found to be well described by a line shape 
of a Gaussian coupled to a Lorentzian, 

y Cfe -f ^l 2 + C^, (9) 



with respect to i'6Si/2:n"i';„ ' Ai"L';„, and Bn"L';„ - Here 
^expt standard uncertainty in f^^^^ determined by 

the fit to the experimental data. For N peaks and all 
three parameters (I'eSi/a:""^';,, ' ^n"L'^„, and Bn"L"j„) 
one would expect the minimum value of the to be 
N — 3. We consistently observed larger values of the 
X^. We attribute this uncertainty to residual misalign- 
ment, which will be discussed in more detail below. In 
order to account for this larger x^ value, we normalize the 
X^ function so that its minimum corresponds to the ex- 
pected value. We then determine the uncertainty in the 
two-photon optical frequency and the hyperfine coupling 
constants by finding the values of i^eSi/2-n"L'^,,, ^n"L'j„j 
and Bn"L" where the normalized x^ function increases 
by one from the point at which it is a minimum. This 
amounts to increasing the uncertainty to make the fre- 
quencies extracted from the different peaks self consis- 
tent. 



B. Results for the 8S1/2 state 

The spectral range of the comb allowed the 8S1/2 state 
to be studied by use of transitions through either the 
6P1/2 (895 nm and 761 nm) or 6P3/2 (852 nm and 795 
nm) states. This redundancy provided a check for con- 
sistency of the measurement and the parameter extrac- 
tion method. A complete study of the systematics was 



where and C/' are the Gaussian and Lorentzian 
amplitudes, respectively, jg and 7^ are the Gaussian 
and Lorentzian widths, fi is the line center of the z*^ 
peak, and is a constant base-line offset. We at- 
tribute the deviation of the experimental line shape from 
the expected Lorentzian line shape to residual misalign- 
ment of the two beams, which we discuss further be- 
low. We designate the resonant repetition-rate frequency 
for the peak determined by fitting the calculated 
spectra for a given i^6Si/2:n"L;;„ , ^wl';,,, and Bn"L';„ 

by /r'''(i'6Si/2:n"L'j„,-4„,/L//^^,B„»L//^^y The resonant 

repetition-rate frequencies determined from the experi- 
mental spectra we designate as f^^P*, The best fit val- 
ues for two-photon frequency and the hyperfine constants 
were determined by minimizing the weighted x^ function 
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FIG. 6: The calculated spectrum (upper trace) and exper- 
imental signal (lower trace) of the 7P — ^ 6S fluorescence 
when the Cs atoms are excited to the 8S1/2 state through the 
6P1/2 state. The peaks in the experimental spectrum that 
are labeled with letters were used in extracting the transition 
frequencies as described in the text. 



performed only for the transitions that occurred through 
the 6P3/2 states. 

Figure [6] shows the spectrum for the 
6S1/2 —>'6Pi/2 ^8Si/2 transition along with the cal- 
culated spectrum. The spectrum was collected by use 
of a 890 nm filter with a 10 nm passband for Fl and 
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a combination of a 755 nm filter with a 40 nm pass 
band with a long-pass filter at 715 nm for F2. The 
13 peaks, labeled a through m in Fig. [6l were used 
to extract the transition frequencies. The peaks were 
acquired with both increasing and decreasing scans of 
/rep, resulting in 26 individual peaks. By use of the pro- 
cedure described above, the extracted center-of-gravity 
two-photon frequency and the hyperfine A constant for 
the 8S1/2 state were found to be 729 009 798.844(38) 
MHz and 219.137(19) MHz. The was 31 for 26 
peaks. Here the uncertainties are determined from the 
function, as described in Section IIV Al and do not 
include any systematic uncertainties. 

The 6S1/2 ^6P3/2 ^^881/2 transition is shown in Fig. SI 
For this spectrum, one of the counter-propagating laser 
beams was filtered with an interference filter centered at 
850 nm with a 10 nm pass band. This resulted in an 
average intensity of ~ 60 W/cm in the interaction re- 
gion, spread over w 4000 optical modes. This light served 
to excite the first step of the two-photon transition at 
852 nm. Based on the transmission profile of the filter 
we estimate that we have « 10 mW/cm^ of light in the 
optical mode resonant with the 6S1/2 ^6P3/2 transition. 
The second stage of the transition was selected by use 
of a combination of a 800 nm short-pass filter and a 
780 nm long-pass filter and provided w 70 W/cm^ over 
w 8000 modes. Using the procedure described above we 
fit the 14 peaks labeled in the figure and find a center-of- 
gravity frequency of 729 009 798.863(29) MHz and a hy- 
perfine coupling constan 219.133(17) MHz, where again, 
the uncertainties are determined purely from the func- 
tion. The value of the was 64 for 28 peaks (14 peaks 
from both the increasing and decreasing frequency scans) 
prior to the normalization. As described in Section [IV Al 
the uncertainties are increased to make the x^ consis- 
tent with the expected value. The resulting frequencies 
are in good agreement with the value extracted from the 
transition through the 6P1/2 state. 

We analyzed in detail the systematic effects that could 
lead to a shift of the measured frequencies from the 
true transition frequencies using the transitions to the 
8S1/2 state through the 6P3/2 state. Systematic shifts 
might be more significant for excitation through the 
6P3/2 state, because of the smaller hyperfine structure 
shifts of the 6P3/2 state, relative to the 6P1/2 state. In 
addition, the signal-to-noise ratio of the spectrum ex- 
cited through the 6P3/2 state was higher than that ex- 
cited through the 6P1/2 state, as an artifact of the filters 
used in the measurements. 

We considered four possible sources of systematic ef- 
fects: ac Stark shifts, Zeeman shifts, pressure shifts, 
and errors arising from misalignment of the counter- 
propagating laser beams. The associated uncertainties 
are summarized in Tab. U and described below. 

To evaluate the systematics we focused on the peaks 
labeled d, e, and / in Fig. S) A typical data set and 
fit to the peaks is shown in Fig. [51 The effect of ac 
Stark shifts on the extracted frequencies was investi- 



gated by varying the optical powers. Neutral density 
filters were used to reduce the power in both beams 
by approximately the same amount. Data were col- 
lected at three optical powers ranging from the maxi- 
mum power, which was used for collection of the full 
spectrum shown in Fig. [31 to 60% of the maximum power. 
The three peaks used in this analysis correspond to 
the 6S1/2 {F = 4) ^6P3/2 {F' = 3, 4, 5) ^881/2 {F" = 4) 
transitions. Because the final state is the same for all 
three peaks, it is not possible to extract both the center- 
of-gravity frequency and the hyperfine A constant from 
these data. In order to place a limit on possible effects 
due to the ac Stark shifts on the center-of-gravity fre- 
quency, we fixed the hyperfine A constant to the value de- 
termined from the analysis of the full spectrum and then 
extracted a center-of-gravity frequency at each power by 
finding the minimum in the x^ function as defined by Eq. 
PH)) . The center-of-gravity frequencies as a function of 
power were fitted to a straight line. We used the slope 
of the line to look for any possible dependence of the 
center-of-gravity frequency on the power. We find the 
slope to be —43(123) kHz/P, where P is the normalized 
operating power. Combining the value of the slope with 
its uncertainty we arrive at a maximum possible shift of 
166 kHz at the the nominal operating power, P = 1. We 
take this value as an estimate of the systematic uncer- 
tainty from the ac Stark effect. Because we fix the hy- 
perfine A constant this approach provides an upper limit 
on the center-of-gravity frequency. Therefore, we do not 
attempt to apply a correction to the data. To determine 
the possible effect on the hyperfine splitting we fix the 
center-of-gravity frequency to the value determined from 
the analysis of the full spectrum and extract the hyper- 
fine A constant from the data. Fitting the data for the A 
constant we found a slope of —22(72) kHz/P, resulting 
in an uncertainty of 94 kHz. 

For two-photon transitions between S states with 
AMp = there is no linear Zeeman shift. However, the 
laser polarization is not perfectly linear and these states 
are step-wise resonant through the intermediate P states, 
which are magnetically sensitive. In addition, the power 
on the first stage of the transition is near saturation, 
so there may be some optical pumping that could couple 
with an external magnetic field to lead to asymmetric line 
shapes. To investigate the effect of Zeeman-shifts due to 
imperfect cancelation of stray magnetic fields, data were 
collected without the compensation coils and also with 
the current through the compensation coils reversed, thus 
doubling the residual magnetic field. Again, the peaks la- 
beled d, e, and / in Fig. [31 were used. The value of the 
hyperfine A constant was fixed and the center-of-gravity 
frequency was extracted as described above. We see no 
evidence for a systematic shift as the magnetic field is 
increased (Fig. [7]). The standard deviation of the data is 
70 kHz, which we take as the uncertainty associated with 
the magnetic field. A similar analysis for the hyperfine A 
constant (with the center-of-gravity frequency fixed) also 
shows no evidence of any shift with the magnetic field. 



8 



Using the standard deviation of the points shown in Fig. 
[7] as an estimate of the uncertainty we assign a 40 kHz 
uncertainty to the A constant due to possible magnetic 
fields. 

For an analysis of the shifts due to collisions, the pres- 
sure of the Cs vapor was changed by heating the vapor 
cell as described above. Data were taken at tempera- 
tures of 297, 318.6, and 345.5 K, corresponding to Cs 
vapor pressures of 0.15, 1.2, and 11 mPa. As described 
above, we determine the center-of-gravity frequency with 
a fixed hyperfine A constant for each pressure. We see no 
indication of any variation of the data with the Cs vapor 
pressure (Fig. [7]). The standard deviation of the frequen- 
cies extracted at the different temperatures is 31 kHz, 
which we take as the uncertainty in the center-of-gravity 
frequency due to temperature or pressure effects. For 
the hyperfine A constant, the scatter in the data taken 
at different temperatures gives a standard deviation of 
18 kHz, which we adopt as the uncertainty. 

The final systematic considered was possible error due 
to imperfect alignment of the counter-propagating laser 
beams. The effect of misalignment of the two beams 
is to broaden and shift the peaks through a first-order 
Doppler shift. As can be seen from Eq. ([2|), the position 
of the peaks depends on the difference in the magnitude 
and direction of the two vectors ki and ^2 as well as 
the velocity class excited. Additional first-order Doppler 
shifts come into play if the wavefronts of ki and ^2 are 
not precisely counter-propagating, which can result from 
divergence of the counter-propagating beams and mis- 
match of the spatial modes. For a spectrum that contains 
many peaks coming from many different velocity classes 
the shift due to the misalignment largely averages out. 
This effect was investigated theoretically by modifying 
the program to include the effect of misalignment of the 
two beams. This was done by changing the /c2 ■ ^ term in 
Eq. ^ to include an additional perpendicular velocity 
component. For the 8S state three spectra were calcu- 
lated with different angles of misalignment. Analysis of 
the calculated spectrum and an estimate of the maximum 
possible misalignment of the beams gives a conservative 
limit on possible effects from angular misalignment of 
22 kHz for the center-of-gravity frequency and 20 kHz 
for the hyperfine A constant. If we attribute the scatter 
in the experimental results to angular misalignment, we 
conclude that for all the spectra, the effective angular de- 
viation due to either misalignment or to divergence and 
mode mismatch of the beams is less than 0.04 radians. 
We note that for the calculated data the misalignment 
of the counter-propagating laser beams leads to no sig- 
nificant shift in the extracted frequency, it does lead to 
a larger spread of the data and an increased value of i 
as observed. This additional spread of the data beyond 
what is expected from the statistical noise is taken into 
account by conservatively increasing the uncertainties to 
achieve the expected x^, as described in Section HV Al 

In addition, the effect of misalignment was investigated 
experimentally on the 6S1/2 —>'6Pi/2 —!>7D3/2 transition. 
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FIG. 7: The dependence of the center-of-gravity frequencies 
for the 6S1/2 — ^6P3/2 — >'8Si/2 state on the power (upper left), 
the magnetic field (upper right), and the Cs vapor pressure 
(bottom). The straight line for the power dependence is a 
weighted linear fit to the data. The vertical axis has been 
offset and centered so that the nominal operating point is at 
zero for the power-dependence plot. The magnetic field and 
pressure data are shown scattered about the weighted mean 
of the data. The horizontal axis for the power dependence is 
plotted as a function of the fraction of the maximum power. 
The horizontal axis for the magnetic field dependence is plot- 
ted as a function of the residual magnetic field, ~ 50 /iT. 



The study of this transition confirmed that while the in- 
dividual peaks do shift, there was no shift in the mean 
value of the transition frequency, within the uncertainty 
of our measurement. 

The non-zero second-order Doppler shift {v/cj'^h'Q is 
about 300 Hz for this transition and is negligible com- 
pared to the other uncertainties. 

Table [J summarizes the corrections and associated un- 
certainties on the extracted center-of-gravity and hyper- 
fine A constants. We see no evidence for any shift that 
is larger than the uncertainty associated with evaluat- 
ing the effect. As our final value we take the mean 
of the 8S1/2 frequencies determined from the 6Pi/2and 
6P3/2 intermediate states, weighted by their respec- 
tive statistical uncertainties, and apply the corrections 
and systematic uncertainties (propagated in quadrature) 
listed in Table HI We find final values for the center-of- 
gravity frequency and hyperfine A constant of 8S1/2 state 
of 

i^6Si/,:8Si/, = 729 009 798.86(19) MHz 
Ass,/, = 219.14(11) MHz. 

These results are in good agreement with previous mea- 
surements of the 8S1/2 state, summarized in Table HIl 



C. Results for the 9Si/2 state 

Figure [5] shows the spectrum obtained for the cascaded 
two-photon transition 6S1/2 — >'6P3/2 — >9Si/2 and the cal- 
culated spectrum. This spectrum was collected by use of 
a 850 nm filter with a 10 nm passband for Fl and a 658 
nm filter with a 10 nm passband for F2. The measure- 
ment of the 9Si/2 state through the 6P1/2 intermediate 
state was not possible because the complementary radia- 
tion at 635.5 nm was on the edge of the laser bandwidth. 

The peaks a-p were used to determine the center-of- 
gravity frequency and the hyperfine A constant as de- 
scribed above. The extracted center-of-gravity frequency 
is 806 761 363.380(58) MHz and the hyperfine A constant 
is 109.932(30) MHz. The was 83 for 32 peaks. 

The systematic effects due to ac-Stark shifts, Zee- 
man shifts, pressure shifts, and misalignment effects 
were investigated in the same way as they were for the 
8S1/2 state. The systematic effects were studied by fo- 
cusing on the three peaks labeled h, i, and j in Fig. HI 
corresponding to transitions from the F = 3 ground state 
hyperfine level to the F" — 3 excited state hyperfine level 
through the F' — 2, 3, and 4 intermediate hyperfine lev- 
els. 

The laser intensity for both transitions was varied. For 
the 6S1/2 ~>6P3/2 stage of the transition the variation was 
from « 40 W/cm^ to 20 W/cm^ and was spread over 
w 4000 optical modes. A linear fit of the center-of-gravity 
frequency to the power gave a dependence of —42(127) 
kHz/P where P is the normalized operating power. We 
therefore assign a 169 kHz uncertainty to the center-of- 
gravity frequency due to possible light shifts (Fig. [S]). 
Similarly, for the hyperfine A constant we found a linear 
dependence consistent with zero slope, 21(58) kHz/P. 
We take 79 kHz uncertainty to be the uncertainty in the 
A constant due to possible ac-Stark shifts. 

As with the 8S1/2 state, no discernable dependence of 
the frequency on the external magnetic field was observed 
(Fig. Based on the scatter (standard deviation) of 
the data we arrive at an uncertainty of 18 kHz for the 
center-of-gravity frequency and 7 kHz for the hyperfine 
A constant. 

Spectra were taken at four different temperatures of 
the vapor cell corresponding to Cs vapor pressures of 
0.15, 1.06, 3.9, and 7.6 mPa. The dependence of the 
center-of-gravity frequency is shown in Fig. [S] There is 
evidence of a possible shift in the frequency with pressure. 
However, even at pressures that are 50 times larger than 
the nominal operating pressure, we see only a 67 kHz 
shift. We again use the standard deviation of the data to 
place a conservative limit on possible uncertainties due 
to pressure shifts and find a 34 kHz uncertainty in the 
optical frequency due to temperature dependent effects. 
A similar analysis for the hyperfine A constant gives an 
uncertainty due to temperature dependent effects of 15 
kHz. 

An analysis of the frequency dependence on misalign- 
ment done in the same way as described above yielded 
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FIG. 8: The calculated spectrum (upper trace) and exper- 
imental signal (lower trace) of the 7P — > 6S fluorescence 
when the Cs atoms are excited to the 9Si/2 state through the 
6P3/2 state. The peaks in the experimental spectrum labeled 
with letters were used in extracting the transition frequencies, 
as described in the text. 
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FIG. 9: The dependence of the center-of-gravity frequencies 
for the 6S1/2 — >'6P3/2 — >'9Si/2 state on the power (upper left), 
the magnetic field (upper right), and the Cs vapor pressure 
(bottom). The straight line for the power dependence is a 
weighted linear fit to the data. The vertical axis has been 
offset and centered so that the nominal operating point is at 
zero for the power-dependence plot. The magnetic field and 
pressure data are shown scattered about the weighted mean 
of the data. The power is plotted as a function of the fraction 
of the maximum power. The magnetic field is plotted as a 
function of the residual magnetic field, « 50 /xT. 



a 58 kHz uncertainty for the center-of-gravity frequency 
and 20 kHz for the hyperfine A constant. 

The uncertainties to the center-of-gravity frequency 
and the hyperfine A constant are summarized in Table HI 

We arrive at a final value for the center-of-gravity fre- 
quency for the 6S1/2 — ^9Si/2 transition and the hyperfine 
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FIG. 10: The calculated spectrum (upper trace) and exper- 
imental signal (lower trace) of the 7P — > 6S fluorescence 
when the Cs atoms are excited to the 7D3/2 state through the 
6P1/2 state. The peaks in the experimental spectrum labeled 
with letters were used in extracting the transition frequencies 
as described in the text. 



A constant of the 9Si/2 state of 



t'6Si/2:9Si/2 
^9Si/2 



806 761 363.38(19) MHz 
109.93(9) MHz. 



This is in good agreement with the previous measure- 
ments of this transition, summarized with other data in 
Table nil and improves on the center-of-gravity frequency 
by two orders of magnitude. 



D. 



Results for the 7D3/2 state 



The center-of-gravity frequency, magnetic dipole, and 
electric quadrupole hyperfine coupling constants for the 
7D3 /2 state were determined from the spectrum in Fig. 
[TUl The filters used in collecting this spectrum were a 
890 nm filter with a 10 nm passband for Fl and a 670 
nm filter with a 30 nm passband for F2. This spectrum 
was taken by exciting Cs through the 6P1/2 intermediate 
state. Use of the 6P1/2 state eliminates possible transi- 
tions to the 7D5/2 state due to angular momentum se- 
lection rules. The groups consist of six peaks, corre- 
sponding to the F -> i^' = 3 ^ F" = 2,3,4 and 
F ^ F' = A ^ F" = 3,4,5 transitions. The 18 peaks 
labeled a through r in Fig. [TU] were used to extract the 
frequency and hyperfine coupling constants. 

The 7D state has J > 1/2 and, consequently, an elec- 
tric quadrupole coupling denoted by the hyperfine B con- 
stant. Minimizing the A, B) function defined by 
Eq. pU]) with respect to the center-of-gravity frequency 
and the hyperfine A and B constants, we find a value 
for the center of gravity of 780 894 762.250(28) MHz, a 
hyperfine A constant value of 7.386(6) MHz, and a hy- 
perfine B constant value of —0.182(50). The was 78 
for 36 peaks. 
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FIG. 11: The dependence of the center-of-gravity frequencies 
for the 6S1/2 — >'6Pi/2 — >-7D3/2 transition on the power (upper 
left), the magnetic field (upper right), and the Cs vapor pres- 
sure (bottom). The straight line for the power dependence is 
a weighted linear fit to the data. The vertical axis has been 
offset and centered so that the nominal operating point is at 
zero for the power-dependence plot. The magnetic field and 
pressure data are shown scattered about the weighted mean 
of the data. The power is plotted as a function of the fraction 
of the maximum power. The magnetic field is plotted as a 
function of the residual magnetic field, ~ 50 /iT. 



The systematic effects were studied by focusing on the 
six peaks labeled g, h, i, j, k, and / in Fig. llOi corre- 
sponding to transitions from the F = 3 ground state. 
The laser intensity on both stages of the transition was 
varied. For the 6S1/2 — ^6Pi/2 stage, the change in inten- 
sity ranged from « 19W/cni^ to « 3W/cm^, spread over 
« 4000 optical modes within the 10 nm (3.8 THz) fifier 
bandwidth. The slope of the center-of-gravity frequency 
as a function of light power yields a dependence of 22(27) 
kHz (Fig. fTTj) . corresponding to an uncertainty of 49 kHz. 
We find uncertainties due to ac-Stark shifts of 14 kHz for 
the A constant and 17 kHz for the B constant. 

No effect of the external magnetic field was observed 
(Fig. [TT|). within the measurement uncertainties. Based 
on the scatter of the data, we assign uncertainties of 17 
kHz, 3 kHz, and 17 kHz for the center-of-gravity fre- 
quency, the hyperfine A constant and the hyperfine B 
constant, respectively. 

Spectra were taken at four different temperatures of 
the vapor cell, corresponding to a Cs vapor pressure of 
0.15, 0.2, 0.3, and 2.0 mPa. As with the 9Si/2 state, there 
appears to be a small dependence of the center-of-gravity 
frequency as a function of the pressure in the cell, lead- 
ing to a 28 kHz shift when the pressure is increased by 
a factor of « 10. However, since we fix the A and B 
constants, this serves as a maximal shift in the center-of- 
gravity frequency and we do not attempt to correct for 
this shift. We take the standard deviation as a conser- 
vative estimate of the error and assign an uncertainty of 
12 kHz to the center-of-gravity frequency due to possi- 



ble temperature-dependent effects. The hyperfine A con- 
stant shows no such dependence and the scatter is less 
than the uncertainty on the mean. We therefore use the 
uncertainty on the weighted mean of the data taken at 
different temperatures to arrive at an uncertainty of 1 
kHz for the hyperfine A constant. For the hyperfine B 
constant we find an If kHz uncertainty (from the scatter 
of the data). 

A theoretical analysis of the frequency dependence on 
misalignment similar to that described for the 8S1/2 state 
analysis yielded a 48 kHz uncertainty for the center-of- 
gravity frequency and uncertainties of 12 kHz for the hy- 
perfine A constant and 143 kHz for the B constant. In 
addition, the effects of misalignment were investigated 
experimentally. Spectra of the six peaks labeled g, h, 
i, J, fc, and I in Fig. [TU] were taken at three different 
beam alignments. The degree of misalignment was char- 
acterized by the peak amplitude. Misalignment of the 
beams so that the peak amplitudes decreased by more 
than a factor of two resulted in no detectable change in 
the center-of-gravity, A constant, or B constant. 

The uncertainties to the center-of-gravity frequency 
and the hyperfine A constant are summarized in Table ID 

With the uncertainties from Table U we find 



i^6Si,2:7D 



J3/2 



780 894 762.250(77) MHz 
7.386(15) MHz 
-0.18(16) MHz. 



These numbers are in good agreement with previous mea- 
surements, shown in Table HI] and lead to a significant 
improvement in the knowledge of the center-of-gravity 
frequency. 



E. 



Results for the 7D5/2 state 



Figure [12] shows the experimental fluorescence spec- 
tra for excitation of the 7D3/2and 7D5/2 states through 
the 6P3/2 state along with the calculated spectra for exci- 
tation to the 7D3/2 and 7D5/2 states. The spectrum was 
collected using a 850 nm filter with a 10 nm passband for 
Fl and a 700 nm filter with a 25 nm passband for F2. The 
close energy spacing of the 7D3/2 and 7D5/2 states did not 
allow for a selection of only one of the states when ex- 
cited through the 6P3/2 state. However, because we were 
able to determine the 7D3/2 state energy and hyperfine 
splitting as described in Section lIVDi we were able to 
fix the frequency of the 7D3/2 state from these measure- 
ments and use the spectrum shown in Fig. [1^ to extract 
information about the 7D5/2 state. In addition, angular 
momentum considerations result in much stronger tran- 
sition probabilities for the 6P3/2 ^7D5/2 transition com- 
pared to the 6P3/2 ~>'7D3/2 transition. To remove the 
effect of the overlap, we calculated the 7D3/2 spectrum 
through the 6P3/2 state using the values determined by 
the 6S1/2 — >6Pi/2 ^7D3/2 excitation. Figure [13] shows 
the contribution of the transitions to the 7D3/2to 
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FIG. 12: The calculated spectrum (upper trace) of the 
7P — >■ 6S fluorescence when the Cs atoms are excited to the 
71)3/2,5/2 state. The experimental signal (lower trace) shows 
the 7P — >■ 6S fluorescence when the Cs atoms are excited to 
the both the 7D3/2, 5/2 states through the 6P3/2 state. In the 
experimental spectrum data are shown for only those regions 
where there was a signal. 



the fluorescence spectrum. We superposed that cal- 
culated spectrum on the experimental spectrum and 
adjusted the peak intensities to match the isolated 
7D3/2 peaks. In this way we accounted for and removed 
the 7D3/2 parameters. To quantify the influence of the 
7D3/2 state and estimate the uncertainty resulting from 
the overlapped spectra, we fit the spectra with and with- 
out subtracting the spectrum from the 7D3/2 state. We 
found a maximum deviation of 10 kHz in the center-of- 
gravity frequency of the 7D5/2 state, 1 kHz in the hyper- 
fine A constant, and 30 kHz in the hyperfine B constant 
which we take to be the uncertainties due to this overlap. 
The contributions from the two transitions are separated 
in Fig. [13] 

For a given ground state hyperfine component F, there 
are nine distinct transitions due to the manifold of in- 
termediate states and upper states. The five groups 
of peaks labeled in Fig. [T^] were used to extract values 
for the center-of-gravity and hyperfine A and B con- 
stants. In group one, only three peaks were of sufficient 
amplitude to be fit. All nine peaks were fit in groups 
three and five. Seven peaks were used in group two 
and six peaks in group four. From the fit to these 34 
peaks we determined the center-of-gravity frequency to 
be 781 522 153.682(25) MHz and the hyperfine A and B 
constants were -1.717(4) MHz and -0.182(86) MHz, re- 
spectively. To improve the precision of the fits, we use 
scans with increasing frequency and decreasing frequency 
simultaneously. The minimum was 90 for 34 peaks 
prior to normalization. 

The systematics were studied in the same way as with 
the previous states. We focused on group four of Fig. 
[T51 due to its size and also the fact that there is very 
little contribution from the 7D3/2 state for this group. 
The optical intensity was varied for both stages of the 
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FIG. 13: The experimental spectrum (upper trace), 
fit to tlie 6S1/2 —^6P3/2 ^7D5/2 transition (middle 
trace), and of the calculated contribution from the 
6S1/2 — >-6P3/2 ^7D3/2 transition (lower trace). Note that the 
vertical scales are different for the two figures. 



transition. For the first stage the intensity over the 10 nm 
bandwidth was varied from w 22 W/cm^ to w 5 W/cm^. 
Data were taken at cell temperatures corresponding to 
saturated vapor pressures of 0.15, 1.05, 3.9, and 8.5 mPa. 
The dependence of the center-of-gravity on the varied 
parameters is shown in Fig. [TJ] and are summarized in 
Table H 

As with the previous transitions there was no signifi- 
cant shift arising from the light power. A linear fit gave 
a slope of —23(62) kHz/P, where P is the normalized 
operating power. This gives us an uncertainty of 85 kHz 
in the center-of-gravity frequency. We found uncertain- 
ties in the hyperfine A and B constants of 11 kHz and 
458 kHz, respectively. 

Unlike the other transitions studied here, we did ob- 
serve a dependence of the center-of-gravity on the mag- 
netic field (Fig. [HI) . There was clear evidence of broad- 
ening in the peaks at the largest magnetic fields. We 
attribute the shift to a break down of the line-shape 
model as the line broadens coupled with the overlap of 
the peaks. In this case, we base our estimate of the uncer- 
tainty on the difference of the point taken at the zeroed 
field with the residual field and find an uncertainty of 
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FIG. 14: The dependence of the center-of-gravity frequencies 
for the 6S1/2 — >'6P3/2 — >-7D5/2 transition on the power (upper 
left), the magnetic field (upper right), and the Cs vapor pres- 
sure (bottom). The straight line in the power dependence 
plot is a linear weighted fit to the data. The vertical axis has 
been offset and centered so that the nominal operating point 
is at zero. The power is plotted as a function of the fraction 
of the maximum power. The magnetic field is plotted as a 
function of the residual magnetic field. 



35 kHz in the optical frequency. There is no evidence for 
a dependence of the extracted hyperfine A constant, and 
we find a magnetic-field uncertainty of 6 kHz for it. For 
the B constant we find an uncertainty of 125 kHz. 

We see no clear evidence for any temperature depen- 
dent effects. Using the larger of standard deviation or the 
uncertainty in the mean, we find uncertainties of 7 kHz, 
1 kHz, and 38 kHz for the temperature-dependent uncer- 
tainties of the center-of-gravity frequency, the hyperfine 
A constant, and the hyperfine B constant, respectively 

Applying these uncertainties we find 

i'6Si/2:7D5/2 = 781 522 153.68(16) MHz 
^7D5/2 = -1.717(15) MHz 
P7D5/2 = -0.18(52) MHz. 



V. CONCLUSIONS 

We have used the technique of velocity-selective two- 
photon excitation to measure the center-of-gravity fre- 
quencies and the hyperfine coupling constants for the 
8S1/2 , 9Si/2 , 7D3/2 , and 7D5/2 states of Cs. The results 
for all of the transitions are summarized and compared 
to other published experimental values in Table HT] This 
technique has potential metrological applications for sec- 
ondary frequency standards at optical wavelengths. At a 
more fundamental level, the hyperfine coupling measure- 
ments are valuable in improving atomic structure calcu- 
lations. The measurement accuracy of multiple transi- 
tions at significantly different wavelengths rivals that of 
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cw laser spectroscopy techniques and demonstrates the 
versatility of direct frequency comb spectroscopy. 

Our measured transition frequencies are in good 
agreement with previous measurements, while achieving 
greater accuracy for the 9Si/2 , 7^3/2 i and 7D5/2 states. 



In these measurements the uncertainties were typically 
limited by signal-to-noise ratio and fitting uncertainties 
combined with uncertainties in the determination of the 
ac Stark shifts and geometrical factors. 



TABLE I: Summary of the contributions to the error budget for the measured transition frequencies and hyperfine constants. 
All frequencies are given in kilohortz. 



Effect 


^'681/2:881/2 


^881/2 


^^681/2:981/2 


^9Si/2 


2^681/2:703/2 


A7D3/2 


-57D3/2 


t'68i/2:7D5/2 


^7D5/2 


-B7D5/2 


ac Stark 


166 


94 


169 


79 


49 


14 


17 


85 


11 


458 


B-Field 


70 


40 


18 


7 


17 


3 


17 


35 


6 


125 


Cs Vapor Pressure 


31 


18 


34 


15 


12 


1 


26 


7 


1 


38 


Alignment 


22 


20 


58 


20 


48 


12 


143 


34 


7 


180 


Overlap of 7D3/2 
















9 


1 


30 


Statistical 


24 


13 


58 


30 


28 


20 


72 


25 


4 


86 


Total 


186 


106 


192 88 


77 


15 


164 


102 


15 


517 



TABLE II: Comparison of the center-of-gravity frequencies and hyperfine coupling constants measured in this work with previous 
measurements. All numbers are in megahertz. The uncertainties represent the standard error (68% confidence interval). 



Measurement 




Center-of-Gravity Frequency 


Hyperfine Constant A 


Hyperfine Constant B 


8S1/2 


This Work 




729 009 798.86(19) 


219.14(11) 




Fendel et al. [22] 




729 009 799.020(26) 


219.125(4) 




Hagel et al. [211 




729 009 798.82(20) 


219.12(1) 




9Si/2 


This Work 




806 761363.38(19) 


109.93(9) 




Weber & Sansonetti 


[23] 


806 761 372(8) 






Farley et al. [29] 






110.1(5) 




Gupta et al. \30] 






109.5(2.0) 




7D3/2 


This Work 




780 894 762.250(77) 


7.386(15) 


-0.18(16) 


Weber & Sansonetti 


[23] 


780 894 772(8) 






Kortyna et al. [311 






7.36(3) 


-0.1(2) 


Auzinsh et al. [32] 






7.4(2) 




7D5/2 


This Work 




781522 153.68(10) 


-1.717(15) 


-0.18(52) 


Weber & Sansonetti 


[23] 


781 522 153(8) 






Auzinsh et al. [3^ 






-1.56(9) 
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